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When scale separation in space and time is poor, the a effect and turbulent diffusivity have to be replaced by integral 
kernels. Earlier work in computing these kernels using the test-field method is now generalized to the case in which both 
spatial and temporal scale separations are poor. The approximate form of the kernel is such that it can be treated in a 
straightforward manner by solving a partial differential equation for the mean electromotive force. The resulting mean- 
field equations are solved for oscillatory a-shear dynamos as well as o? dynamos in which a is antisymmetric about the 
equator, making this dynamo also oscillatory. In both cases, the critical values of the dynamo number is lowered by the 
fact that the dynamo is oscillatory. 
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1 Introduction 

Mean-field dynamo theory describes the evolution of the av- 
eraged magnetic field. This theory is relevant for the under- 
standing of the origin of ordered magnetic fields in the Sun 
and other late-type stars. Compared to the original induction 
equation, the averaged equation contains extra terms which 
capture the effects of systematic correlations between veloc- 
ity and magnetic field fluctuations. Some of these terms (for 
example the a effect) can be responsible for the generation 
of mean magnetic fields. 

Mean-field dynamo theory provides an important tool 
for a number of astrophysical applications. However, it also 
suffers from several shortcomings, some of which can be 
the result of simplifications that are not well justified and 
often not even necessary. In this paper we focus on the 
issue of poor scale separation in space and time. Broadly 
speaking, if there is poor scale separation, multiplications 
with mean-field coefficients must be replaced by convolu- 
tions with corresponding integral kernels. Obviously, as far 
as temporal scale separation is concerned, this effect can- 
not be very important for the Sun, because the cycle time 
is much longer than the convective turnover time. However, 
with respect to spatial scale separation this is no longer true, 
because at the bottom of the solar convection zone the pres- 
sure scale height and with it the typical size of the convec- 
tion cells is 50 Mm, and hence comparable to the depth of 
the convection zone of 200 Mm which is also the scale of 
the mean magnetic field. Although the concept of writing 
the mean electromotive force as a spatio-temporal convo- 
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lution with the mean magnetic field was well known (e.g., 
Radler 1976), there was the problem that, until recently, not 
much was known about the form of the integral kernels that 
are to be used. In the past there have been several attempts 
to compute the integral kernels from turbulence simulations 
(e.g., Miesch et al. 2000; Brandenburg & Sokoloff 2002), 
but the situation has changed drastically with the advent of 
the test-field method (Schrinner et al. 2005, 2007) which 
allowed an accurate determination of the integral kernels 
in space (Brandenburg et al. 2008) and time (Hubbard & 
Brandenburg 2009). As a result, we now know that the ker- 
nels of most of the components of the ol and r] tensors are 
Lorentzians in spectral space and exponentials in real space 
(Brandenburg et al. 2008). It turns out that in these sim- 
ple cases, the resulting integro-differential equation for the 
magnetic field can be reformulated into a set of two coupled 
differential equations of parabolic type, one for the mag- 
netic field and one for the electromotive force. 

In hindsight, we can say that even temporal scale sepa- 
ration can sometimes be relevant, because nowadays we are 
not only comparing with the Sun and other astrophysical 
bodies, but also with direct numerical simulations (DNS). 
In DNS we may well have situations in which the dynamo 
e-folding times and perhaps also the cycle periods become 
comparable to the turnover time of the turbulence. In these 
more extreme situations, we have much better possibili- 
ties of testing theory. Furthermore, with DNS there is more 
freedom in constructing cases that may be hard to find in 
real astrophysical bodies, but for which the same mean- 
field theory should equally well be applicable. Furthermore, 
DNS allow us to determine turbulent transport coefficients 
to high accuracy, facilitating therefore detailed comparison 
with mean-field theory. Indeed, it turns out that in DNS 
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the growth rates of dynamos can well be comparable to 
the turnover time. A dramatic example was presented by 
Hubbard & Brandenburg (2009), where the growth rate of a 
Roberts flow dynamo was found to be significantly different 
from the value expected based on the analytic dispersion re- 
lation using coefficients that have been determined numeri- 
cally using the test-field method, but under the assumption 
of perfect scale separation in time. 

2 Formalism 

To set the scene, let us begin with the mean-field dynamo 
equation for the mean magnetic field B, 



dB 
'dt 



Vx{UxB + £- ?//io J) , 



(1) 



where U is the mean velocity, £ is the mean electromotive 
force, and J — "V x B //io is the mean current density, with 
fiQ being the vacuum permeability, and 77 the microscopic 
(molecular) magnetic diffusivity. Under certain conditions, 
£ can be expanded in terms of the mean magnetic field and 
its derivatives as 



(2) 



where the comma denotes partial differentiation and the 
dots refer to higher spatial derivatives of B, temporal 
derivatives of B, as well as terms independent of B. 

In many cases of practical interest, only the lowest (in- 
cluding the zeroth) order spatial derivatives are retained, be- 
cause they are sufficient for capturing qualitatively new ef- 
fects such as large-scale dynamo action. This has led to a 
large number of mean-field dynamo models that were ap- 
plied to the Sun, other stars, galaxies, and even accretion 
discs, hi such models, the length scales of the resulting 
mean field become often quite small, especially in the non- 
linear regime; see, e.g., Chatterjee et al. (2011a, Figs. 9- 
1 1). In this context, 'small' means that the scale of the mean 
field becomes comparable to and even smaller than the scale 
of the energy-carrying eddies. In stratified turbulence, as 
present in the Sun, the scale of these eddies is often as- 
sumed to be proportional to the local pressure scale height, 
which is about 50 Mm at the bottom of the solar convection 
zone. However, in mean-field models the magnetic fields 
show frequently variations on scales much smaller than this. 
Chatterjee et al. (2011a), discussed the small-scale fields at 
the bottom of the convection zone in their simulations of 
a mean-field dynamo model as an artifact of the neglect of 
nonlocality in space, but no solution to this problem was 
feasible at the time. 

Looking at Eq. (|2]i, it is clear that higher spatial deriva- 
tives need to be retained when the mean field is no longer 
slowly varying in space. Unfortunately, such a series expan- 
sion becomes easily quite cumbersome, and it is then better 
to replace Eq. (|2]i by a convolution of the mean magnetic 
field B with some integral kernel. As alluded to above, a 
representation of £ in terms of a convolution of B with a 



kernel determined by the statistical properties of the turbu- 
lence has long been known to be the more basic one (e.g., 
Radler 1976). By allowing the convolution to be also over 
time, we can automatically include all temporal derivatives 
as well, i.e., we can instead of Eq. (|2|i write 



(3) 



£,{x,t)^ j Gij{x,x',t,t')Bj{x',t')d^x' dt', 

where we have again ignored terms that are independent of 
B. 

For simplicity, we shall restrict ourselves now to statis- 
tically homogeneous and steady turbulence, in which case 
Gij is translation invariant in space and time and depends 
thus only on the arguments x — x' and t — t'. In cases with 
boundaries, this is not possible, but the formalism presented 
below can easily be adapted to such cases as well; see Chat- 
terjee et al. (2011b). 

Continuing now with the translation invariant case, the 
convolution becomes a multiplication in Fourier space, i.e., 

£,ik,u)^G,,{k,Lu)Bj{k,u), (4) 

where hats indicate Fourier transformation in space and 
time, e.g.. 



(5) 



£^ik,u}) = J £,(a;,Oe'^'=-^-"*) d^xdt. 

In view of the traditional distinction of contributions to £ 
from the a effect and turbulent diffusivity, it is convenient 
to write the Fourier transform of the kernel in the form 



Gij{k,uj) 



D{k,uj) 



(6) 



where a\^^ and tjI^I are assumed to be tensors that are inde- 
pendent of k and ut. The goal of this paper is to verify the 
approximate validity of Eq. (|6]l and to consider the conse- 
quences of such a structure for mean-field dynamo models. 

In the following we consider triply periodic domains 
and define mean fields as planar averages over the x and 
y directions, so that B is only a function of z and t. In that 
case, k = (0, 0, fc) has only one component. Recent work 
of Hubbard & Brandenburg (2009) has already revealed that 
for fixed k, D{k, lu) is proportional to 1 — iwr, where t is 
a fit parameter that is approximately equal to the turnover 
time, i.e., ruims^f ~ 1. where u^„is is the rms velocity 
of the turbulence and kf is the wavenumber of its energy- 
carrying eddies On the other hand, for a; = 0, D{k,uj) 
is approximately proportional to 1 + (ak/kf)'^, where a is 
a dimensionless parameter, for which values between 0.2 
and 1 have been found over a range of different simulations 
(Brandenburg et al. 2008, 2009; Madarassy & Brandenburg 
2010). Consequently, we propose in the present paper that 
D{k,uj) can be approximated by 

D{k,Lj) = l-iujT + e^k'^. (7) 

with the additional parameter £ having the dimension of a 
length. Such a form, even if it is still an approximation that 
neglects higher powers of k and oj, has aheady the advan- 
tage of alleviating problems of unrealistic variations of the 
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magnetic field on short length and time scales. Moreover, it 
leads to an easily treatable partial differential equation for 
S in real space, namely 



1 



d 



dt 



dz 



Note that in the limit r — ^ and — > 0, the usual dynamo 
equations are recovered. Thus, nonlocality is captured sim- 



ply by specifying t and £, while the tensors a. 



(0) 



and v-fk 



can simply be regarded as the usual ones of a effect and 
turbulent diffusivity for the limit k ^ 0, uj ^ 0. Therefore, 
the superscripts (0) will from now on be dropped. The pur- 
pose of this paper is to establish not only the validity of this 
approach, but also to assess the properties of mean-field dy- 
namos when £ is obtained as the solution of the evolution 
equation Eq. 

3 The kernel function D{k,u) from DNS 

3.1 Turbulence in a periodic domain 

In the following we present results for three-dimensional 
isothermal turbulence that is being forced in a narrow range 
of wavenumbers around a representative wavenumber k[. 
We adopt a cubic domain of size L^, measure length in units 
of the inverse minimal box wavenumber ki — 2tt/L and 
choose ki/ki « 2.2. We vary the magnetic Reynolds num- 
ber. 



-Km = Ur^s/rjkf, 



(9) 



where Uj-ms is the rms velocity of the turbulence, keeping 
the rms Mach number, Urms/cs at around 0.1. In agreement 
with the considerations above, time is expressed in units of 
the turnover time, defined here as ro — (urmsfct)^^, and 
the turbulent magnetic diffusivity is expressed in units of 
~ Uyj-ns/3kf (cf. Sur et al. 2008). 

3.2 Test-fields in space and time 

To establish the form of Eq. (|7]) we use the test-field method, 
i.e., we solve, for a given turbulent velocity field, the equa- 
tions governing the departure of the magnetic field from a 
given mean field, that is, we determine the magnetic ^McfM- 
ations b caused by the interaction of the turbulent velocity 
with the mean field. This mean field is referred to as the test 
field and is marked by the superscript T. For each test field 
B , we find a corresponding departure 6^ = V x by 
solving the inhomogeneous equation for the corresponding 
vector potential a^, 



I h II I! II h 



dt 



(10) 



where (^u x ti^Y = u x — u x b"^ is the fluctuating 
part of M X 6^, and compute the corresponding mean elec- 
tromotive force, £ = u x b^. We use test fields that are 
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Fig. 1 b{k, uS) for wro = 1.04 (Run A), 0.52 (Run B), 
and 0.26 (Run C). Open and filled circles denote the real and 
imaginary parts of D{k,uj) as obtained from the test-field 
method; the parabolas give a fit proportional to 1 + £^k^. 
Dashed lines: average of the three data points of the imag- 
inary part of D{k,uj) for each Rm- For the fit parameters £ 
and T see Table [T] 



harmonic functions with wavenumber k and frequency u 
and point either in the x or in the y direction, i.e., 

B^'^'^'^— Bi cos kz cos uit, B^^'^'^— Bi sink z cos ojt, (11) 

i = 1,2, where ei and 62 are unit vectors pointing in the x 
and y directions, respectively. The third component is here 
without interest, because V ■ B = dBz/dz — 0, so = 
const, and is chosen to be zero initially. 

Using the standard test-field method, we obtain directly 
the tensors aij{k, uj) and %jfe(^, From that we can de- 
termine D for different values of k and oj according to 

D{k,uj) = a,j{0,Q)/a,j{k,uj), (12) 

or 

D{k,uj) = fj,j{0,0)/fi,j{k,oj), (13) 

employing the known values ctij (0, 0) or f/ij (0, 0). Further- 
more, since we consider isotropic turbulence, both tensors 
are isotropic, i.e., ctij — aSij and T]ijk — T]t£ijk, but a = 
for non-helical turbulence (Runs A-D of Table [T]). For this 



Table 1 Summary of fit parameters for runs without he- 
licity (Runs A-D) and one with helicity (Run E) using either 
Eq. O or Eq. O- 



Run 


Rm 


UJTO 


r/ro 


£ki 


Equation 


A 


8 


1.04 


1.85 


0.99 




B 


8 


0.52 


1.46 


0.88 




C 


8 


0.26 


1.24 


0.83 


{H 


D 


53 


0.38 


1.21 


0.77 


Gl 


E 


57 


0.35 


0.67 


0.60 




E 


57 


0.35 


0.59 


0.80 


03 
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case D{k,uj) is shown in Fig.[Tl where we plot its real and 
imaginary parts for the scale separation ratio fcf/fci = 2.2, 
Rm = 8, and three values of ujtq. The real part of D{k, w) 
is a fit to a profile of the form 1 + ^^fc^, while the imaginary 
part of D{k,uj) is approximately independent of k. This is 
consistent with 3{Z?(fc,aj)} ~ —lot, where r is obtained 
by taking the average value of lot for all three k values. We 
find that t/tq and ikf are of the order of unity. In agreement 
with the ansatz O the parameter £fcf varies only weakly 
with Lo, but t/tq shows a stronger variance, indicating the 
presence of higher powers of lo in D. Both parameters vary 
somewhat with i?^; see Table [T] for details. The additional 
Run E differs from Run D only in including helicity in the 
forcing and hence in the flow. For both t and t the resulting 
values obtained by using Eq. (fT2l l and Eq. (fT3] l are similar 
The value of I is also similar to that of Run D whereas t is 
reduced by a factor of w 2. 

Comparing the results for Runs D and E suggests that 
in Eq. (|7]l the values of r are reduced by a factor of 2 when 
there is helicity in the turbulence, while t remains approxi- 
mately unchanged. 

Thus, in conclusion, we have verified that, for a turbu- 
lent flow such as that considered here, the integral kernel 
in Eq. (|6]l with D{k,uj) is roughly given by Eq. (|7]l. We 
concede, however, that the modeling of the lo dependence 
of Gij is worth to be improved taking into account higher 
orders in oj. In the remainder of this paper we examine prop- 
erties of the resulting mean-field equations. 

4 Application to mean-field dynamo models 

4.1 Nonlocality in dynamo waves 

Some limited insight into the effects of nonlocality for dy- 
namo waves has already been provided in the paper by 
Brandenburg et al. (2008), who considered nonlocality in 
space, but not in time. Based on their test-field results, they 
found a kernel compatible with a Lorentzian in k space. 
Generally speaking, such a kernel makes the resulting mean 
electromotive force smoother by acting preferentially on the 
largest scale in the domain. In the present paper we repeat a 
similar experiment, but with the difference that we include 
here also nonlocality in time. 

Nonlocality in time can lead to somewhat unexpected 
behavior of oscillatory dynamos of a-shear type in that it 
enhances their growth rate and, more importantly, it lowers 
the critical value for dynamo action. This is different from 
the dynamo case, where the presence of an extra time 
derivative always leads to a lower growth rate (Brandenburg 
et al. 2008). This can be seen by comparing the two disper- 
sion relations for and a-shear dynamos with constant a 
and shear By making an ansatz of the form 

B = Bexp [i{kz - cot) + At] , (14) 

with real coefficients k (wavenumber), lo (frequency), and 
A (growth rate), we can easily obtain the dispersion rela- 
tion for the system of Eqs. ([T]) and ^ in implicit form. In 
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Fig. 2 Critical dynamo number for an a-shear dynamo as 
a function of rrytfc^ for different values of £k. Microscopic 
magnetic diffusivity, rj, is here neglected. 

the case of an a^ dynamo with 77 = we have (see Ap- 
pendix |A]T]) 

X = {±\ak\ - ritk^) , w = 0, (15) 
where we have introduced the correction factor 

^^l + TX + fk"^, (16) 

In the case of a pure a-shear dynamo, with the aB^ term 
neglected in favor of SB^, and again ry = 0, it is convenient 
to seek marginally excited oscillatory solutions with A = 0, 
which gives, with ^ = 1 + 



2 ' 



which allows us to compute the critical dynamo number 

D„it = aS/rf^k^ = 2w(l - lo'^t /i^tk") htk" . (18) 

Note that for an a^ dynamo the threshold remains un- 
changed, while for an a-shear dynamo the term lo'^t always 
lowers the threshold. 

In Fig. |2] we show -Dciit for an a-shear dynamo as a 
function of Trjtk^ for different parameters £fc. For small val- 
ues of r, the usual value of -Dcrit = 2 is obtained (e.g., 
Brandenburg & Subramanian 2005). 

4.2 Nonlocality and boundaries 

We have mentioned in the beginning that the effect of spa- 
tial nonlocality should consist in a spatial smoothing of 
the mean electromotive force. However, the solutions pre- 
sented so far are all entirely harmonic. To see the anticipated 
smoothing effect, we can either consider nonlinear solutions 
(as done in Brandenburg et al. 2008), or we can consider so- 
lutions with boundaries, which breaks the monochromatic 
nature of the solutions. 

In the following we solve Eqs. ([U and ^ numeri- 
cally in terms of the mean magnetic vector potential A, so 
B = V X A We use here the PENCIL COD^H which is 
a high-order public domain code (sixth order in space and 



|http : /7pencil-code ■ googlecode . com/| 
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third order in time) for solving partial differential equations, 
including a range of different mean-field equations. The fi- 
nal set of equations is for vanishing mean flow 



dt 



de 

'dt 



aB — rjtJ 



(19) 



(20) 



of which only the x and y components are relevant. We have 
introduced here the additional parameter ry^ = l"^ /t having 
the dimension of diffusivity. 

For a simple dynamo with boundaries we choose the 
dynamo with a linear a profile, 

a(z) = aoz/L^, (21) 

with < z < Lz, where = 7r/2fci is the size of the do- 
main. For the sake of simplicity we retain here the assump- 
tion of isotropy, although it is strictly not tenable under in- 
homogeneous conditions. Here, fci is the lowest wavenum- 
ber for a quarter-cosine wave obeying the boundary condi- 
tions 



Ax.z — An,, — Sx.z — ^-u.z — on z — 0, 



and 

Ax — Ay — Sx — — 



on z = Lz. 



(22) 



(23) 



These conditions correspond to a perfect conductor condi- 
tion on z = and select solutions B antisymmetric about 
z = 0. 

We recall that the dynamos with the linearly varying 
a profile ( |2TI ) are always oscillatory with dynamo waves. 
This was first noticed in direct numerical simulations (Mi- 
tra et al. 2010), but was then also confirmed for mean-field 
models (Brandenburg et al. 2009) and is consistent with the 
parametric survey of solutions given by Riidiger & Holler- 
bach (2004). 

We have computed marginally excited dynamo solu- 
tions for different values of rrj^ kf and r/^ / rjt . For compari- 
son, for the value kf/ki = 2.2 considered in Sect. 13.11 we 
have, using r/t = Vto ™d assuming t/tq = £k{ = 1, 



T77tfc? = {3kf/kl] 



0.06, 



(24) 



The critical values of the dynamo number Cq — a/rjtki and 
the resulting normalized cycle frequencies uj /rj^kf are given 
in Table |2] For five particular cases, denoted by the labels 
(a)-(e), the corresponding butterfly diagrams are shown in 
Fig. El 

Similar to the a-shear dynamos discussed in Sect. 14.11 
we find that the critical dynamo number C"'*, is lowered 
in all cases with r 7^ 0; see Table |2] Furthermore, and 
perhaps somewhat surprisingly, we find that, as rj^/rit is 
increased, the dynamo wave weakens significantly before 
reaching the equator; see panels (b)-(e). On the other hand, 
increasing rrjtkf from 10"'^ to 1 does not affect the weak- 
ening of the dynamo wave near the equator, but it rather 
enhances its speed. Whether similar results also apply to 



Table 2 Dependence of C"'' and normalized cycle fre- 
quency Lj/rjtkl on Trjtkl and rje/rjt for marginally excited 
solutions of dynamos with linear a profile (|2T]) . 



Run 


rritkl 




^crit 


oj/ritkf 


(a) 


0.001 


0.001 


5.16 


1.64 




0.1 


0.001 


4.65 


0.74 


(b) 


1 


0.001 


2.76 


0.88 




1 


0.1 


1.11 


0.87 


(c) 


1 


0.3 


2.84 


0.86 




1 


0.7 


3.68 


0.78 


(d) 


1 


1 


5.30 


0.64 


(e) 


0.06 


3 


8.12 


0.58 



a-shear dynamos is however not obvious. Also, while the 
anticipated smoothing effect might explain the weakening 
of the dynamo wave near the equator, it does not seem to 
operate in the same way in the proximity of the boundary 
at z — Lz- Instead, we see that the dynamo wave is now 
more nearly perpendicular to that boundary compared with 
the case rj^ Q. 

5 Conclusions 

The present work has established that the Fourier transform 
of the integral kernel for the representation of the mean elec- 
ti'omotive force in the isoti^opic case is well approximated 
by 

G{k,uj)G^- : (25) 

which, in turn, can be captured by solving a partial differ- 
ential equation for the mean electromotive force with a first 
order time derivative and a Laplacian that plays the role of a 
diffusion term. Our work has illustrated the great ease with 
which nonlocality in space and time can be implemented in 
a dynamo model. Indeed, the chosen, simplest possible ker- 
nel leads to a rather plausible representation of the partial 
differential equation governing the evolution of the elec- 
tromotive force. Furthermore, the application to spherical 
and other coordinate systems is quite straightforward and 
already fully functional in the PENCIL Code. 

It turns out that, while nonlocality normally hampers 
dynamo action, it can actually make the dynamo more ea- 
sily excitable provided it is oscillatory. This has here been 
shown analytically for standard dynamo waves in the pres- 
ence of shear, but it has also been found in the case of an o? 
dynamo in spherical geometry where the oscillatory behav- 
ior is a consequence of the spatial antisymmetry of a about 
the equatorial plane (Mitra et al. 2010). 

Another issue that has not been addressed here is the 
question of nonlinearity. Our present approach is easily ex- 
tendable to the case where a and r\t are nonlinear functions 
of B, as in the case of usual algebraic quenching. Even 
the case of dynamic a quenching (Kleeorin & Ruzmaikin 
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Fig. 3 Butterfly or zt diagram of By for mean-field mod- 
els with different combinations of Tr]tk\ and rielrix.. (a): 
Tr]tk\ = r]e/rit = lO'^; (b) - (d): TTjtkf = 1, Ty^/r^t = 
10-3, 0.3, 1, (e): rijtkl = 0.06, ?7^/?7t = 3. Tcyc - cycle 
period. 



1982) could easily be included. Here, yet another differ- 
ential equation is being solved, namely one for a magnetic 
contribution to a. One might imagine that the effects of this 
additional equation are already captured by the d£ / dt equa- 
tion. However, it should be remembered that the dynamic a 
quenching also contains effects of magnetic helicity fluxes 



and is capable of reproducing the resistively slow saturation 
in the absence of such fluxes. 

We regard the approach of solving a partial differential 
equation for 5 as a natural one, which supersedes the usual 
dynamo equations where t — j- and ^ — ;> is assumed. 
In many typical situations, neither of the two assumptions 
are well satisfied. We also recall that the approach of in- 
cluding the time derivative of S addresses the problem of 
causality, i.e., the propagation speed of disturbances of B is 
automatically limited to the value of the rms velocity of the 
turbulence, as demonstrated in Brandenburg et al. (2004). 
Furthermore, the presence of the diffusion operator in the 
evolution equation for £ is natural and advantageous be- 
cause it ensures numerical stability and, more importantly, 
it prevents, in a physical way, the emergence of artificially 
sharp structures on scales comparable to or below that of the 
turbulence. Is should be noted, however, that, while the time 
derivative of £ emerges as a natural consequence from the 
T approach (Blackman & Field 2002, 2003), there does not 
seem to be a likewise natural motivation for the presence of 
the diffusion term in the equation for £. 
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A Dispersion relations for nonlocal dynamos 

A.l c? dynamos 

We begin by writing the governing equations (|Tll-l[8j in component 
form for homogeneous turbulence, i.e. constant mean-field coeffi- 
cients, hence 

at az az-^ 

f^^+f^+.^, (A2) 
at az oz'' 

— d£x fid'^S-x , dBy 

£y + r^-e'^=^By-.,t§^. (A4) 
ot oz^ oz 

The dispersion relation is easily obtained by employing the ansatz 

il4i with IX) = in these equations and writing them in matrix 

form, Mq = 0, where q — (Bx, By, Ex, £y)'^ is the state vector 

and 

(X + rik^ +iA; 

\ + rjk^ -ifc 

-Q -ikr]t l + Xr + fk^ 
V +ikrit -a 1 + Ar + fk'^ ' 

is the matrix M for the dynamo. Nontrivial solutions have van- 
ishing determinant, which yields 

[{X + rjk^) (l + Ar + ^'^fc^) +Vtk^Y = a^'^^- (A5) 
Taking the square root leads to the implicit solution il5i . 



A.2 a-shear dynamos 

In the case of a pure Q-shear dynamo with a mean flow of the form 
U — (0, Sx, 0), and neglecting the terms aBx, we have 

dBx dly , d^Bx 
Ot Oz Oz^ 

dBy _ .-^ ,d£x ,d^By 

— — — >bi(a; + — h ?7 „ „ , (A/) 

at oz oz'' 

d£x fid'^Ex , dBy 
£x+T—-~i——-^+'qt——, (A8) 
at oz-' oz 

, d£y ^2d^£y ^ dBx 
£y + T—f--l—-f=aBy-r)t-^. (A9) 
at az-' oz 

In the marginally excited oscillatory case, the matrix M is for 77 = 





M 













(Bx 


-S 




-ifc 




By 





-ikr)t 




:•) 




ikr]t 


— Q 










(AlO) 



^ = 1 + i k . The dispersion relation becomes 

[-M^ - iuJT) + ritk^] ^ + ikaS{i - iuT) = 0. (All) 
Solving separately for real and imaginary parts, we obtain 

- ijj'^i^ + (-cjV + ritk'^f + ■qlk'^ + kaSuiT = (A12) 
and 

2aj(a;V-r7tfc^) + fcaS = 0. (A13) 
Eliminating kaS yields then Eqs. l ll7t and l llSt . 
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